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Abstract 

We compute the real and imaginary parts of the one-loop electroweak contributions to the 
left and right tensorial anomalous couplings of the tbW vertex in the Standard Model (SM). 
For both tensorial couplings we find that the real part of the electroweak SM correction is 
close to 10% of the leading contribution given by the QCD gluon exchange. We also find 
that the electroweak real and imaginary parts for the anomalous right coupling are almost 
of the same order of magnitude. The one loop SM prediction for the real part of the left 
coupling is close to the 3o discovery limit derived from b sj. Besides, taking into account 
that the predictions of new physics interactions are also at the level of a few percents when 
compared with the one loop QCD gluon exchange, these electroweak corrections should be 
taken into account in order to disentangle new physics effects from the standard ones. These 
anomalous tensorial couplings of the top quark will be investigated at the LHC in the near 
future where sensitivity to these contributions may be achieved. 
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1 Introduction 



Top quark physics at the Large Hadron Collider (LHC) is an important scenario for testing 
physics above the electroweak scale. No deviation from the predictions of the Standard 
Model (SM) is found in top data III, nowadays dominated by the Tevatron experiments. 
This situation may change with the LHC already running and taking data. The SM dominant 
decay mode t — t- bW^ will be precisely measured at the LHC and sensitivity beyond tree level 
SM will be achieved. It is generally believed that, due to its large mass, physics of the top 
quark will be useful to probe new theories above the electroweak scale |l2l|3l. At the LHC, 
top new physics may show up in new top quark decay channels or in the measurement of 
the top standard and anomalous couplings [|3l|4l. In renormalizable theories, the anomalous 
couplings appear as quantum corrections, as it is the case for the SM, and also in many new 
physics theories. In a model independent approach there are two ways of parameterize the 
unknown physics at high scales. One is the effective Lagrangian method [5| which is a way 
to describe low energy physics effects originated at a higher energy scale. These effects 
are parameterized with non-renormalized terms invariant under the SM gauge symmetry 
SU (3)c X SU {2)l xU{1)y and written in terms of the low energy (standard) particle spectrum 
fields. It is assumed that the new particles spectrum lies at an energy scale well above 
the electroweak scale. The other way is just by assuming the most general form of the 
Lorentz structure for the tbW amplitude. There are many terms in the effective Lagrangian 
that may give contributions to the same Lorentz structure in the vertex, in particular to the 
tensorial couplings we are interested in. Besides, some of those terms can be rewritten by 
using the equations of motion, so the identification of the effective Lagrangian terms with 
the form factors is not direct nor unique. In this paper we will use the second approach, 
parameterizing in the most general way the tbW amplitude. This approach has the advantage, 
over the effective Lagrangian approach, that it does not break down even if any relatively 
light particles, as new scalars, for example, come into the game. 

Some effects related to the top anomalous couplings, both in the t — )■ bW^ polarized branch- 
ing fractions -for the three helicity W possibilities- and in single top production at the Teva- 
tron and at the LHC, have already been studied in the recent years [|3]|6l. However, at the 
LHC it will be possible to have new suitable observables in order to perform precise mea- 
surements of the anomalous couplings. Some aspect of this top quark physics have also been 
investigated in models with an extended Higgs sector, technicolor models, supersymmetry 
models and Little Higgs models fTl . 

The anomalous couplings are gauge invariant quantities so one can think of testing the SM 
predictions and new theories through observables that are directly sensitive to them. In fact. 
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not only top branching fractions and cross sections are predictions of the models to be con- 
fronted with data but, in the same way as in the past the anomalous magnetic moment for the 
electron gave the first success of quantum field theory, also these tbW gauge invariant tenso- 
rial couplings can be used to check the predictions of new physics theories. One loop QCD 
and electroweak contributions to the top branching fractions for polarized W's have been 
studied in the frame of the SM flSl. These contributions, and the corresponding measure- 
ments, have no special sensitivity to the anomalous couplings which enter in the observables 
as small corrections. The explicit dependence of the polarized branching fractions on the 
anomalous couplings have been computed in [|9lll0||, where also the sensitivity to them has 
been considered. Specific observables that are directly proportional to the tensorial couplings 
have been studied in |[6l [T0l - [T2l and more recently new observables were presented in llT3l . 

In this paper we compute the electroweak SM contribution to the left and right "magnetic" 
tensorial couplings of the tbW vertex and we discuss their observable effects; we also com- 
pare this contributions with some new physics predictions considered in the literature. These 
CP-conserving pieces of the tbW vertex are different from zero only at one loop in the SM 
and the same is true in many extended models. The QCD gluon-exchange contribution to the 
tensorial couplings is the dominant one and has been reported in the literature only for the 
right coupling [I141I15I . The left tensorial coupling is proportional to the bottom quark-mass 
due to the chirality flipping property of this coupling and by the fact that it only couples to 
a right b-quark. For these reasons it is suppressed and it is generally assumed to be negligi- 
ble. However, the measurement of both of them appears as feasible in dedicated observables 
computed for top production at the LHC. 

For the right tensorial tbW coupling the comparison with the SM is usually performed in the 
literature by taking as a reference only the one loop QCD contribution. The most promising 
new physics models predict a few percents deviation from this QCD-value. However, as 
it is shown in this paper, we found that the electroweak contribution is also at the level of 
10% with respect to the leading gluon exchange, and should be taken into account when 
comparing with data. Detailed studies will be necessary in order to disentangle new physics 
contributions from the electroweak standard ones. 

In section 2 we define the anomalous couplings and we review the theoretical as well as the 
experimental status for the physics for which they are involved. In section 3 our computation 
is presented and in the final section we present our conclusions. 
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2 The tensorial tbW vertex: experiment, SM and beyond 



For on-shell particles, the most general amplitude for the decay t{p) — b{k) W^{q) 
can be written in the following way: 



Hbw+ = — — — T^e''*"^ 



{VlPl + VrPr) + '-^^ {elPl + ErPr) 



Ut, (1) 



with Pi^R = (l^js)/!; p,k and q = p — k denote the top, bottom and W boson four momenta, 
respectively. The tensorial left and right magnetic moments are gi and gR respectively. The 
tree level SM couplings are Vr = 0, = Vtb (the Cabibbo-Kobayashi-Maskawa, CKM, 
matrix element), gR = and gi = 0. This expression for the amplitude, written in terms 
of the most general form factors, is appropriate for a model independent analysis of the tbW 
amplitude. The anomalous form factors gR and gi are chirality flipping and dimensionless 
functions of q^. For all particles on-shell, as can be assumed for the top decay, we have 
q^ = M^. Besides, these dipole moments are gauge independent quantities and may be 
measured with appropriate observables. 

These form factors are generated by quantum corrections in the SM. In renormalizable the- 
ories, such as some extensions of the SM, Vr can appear at tree-level while tensor couplings 
gR and gi, are induced as one loop quantum corrections. Values of \Vi\ different from the 
ones given by the global fit [[T6l in the SM Vtb — 1^ that we assume, are not experimentally 
excluded ifTTl and they are still an open window to test new physics. This issue (and also 
possible deviations of Vr — 0) will not be the object of our work, where we concentrate only 
on gR and g^. 

In renormalizable theories, the tensorial couplings are finite quantum corrections quantities 
that do not receive contributions from renormalization counter-terms at one loop. In addition, 
contrary to what happens for the Vr form factor, the tensors gi R couplings are infrared safe 
quantities. One loop QCD corrections generate the leading contribution to the tensorial 
couplings gR and g^. This one loop QCD gluon exchange contribution to gR was computed 
in ifTSl and they found the value g^'^^ = —6.61 x 10^^. Direct observables with sensitivity 
to gR will be accessible to the LHC experiments as was discussed in IfTTl . The left tensorial 
coupling term couples a right b-quark and thus it is proportional to m^,. For this reason, it 
is generally believed that the SM value for gi is much smaller than the one for gR. We will 
show that this in not exactly the case. 

New physics signals can also show up in the analysis of the top decay t — )■ bW^. In particular, 
significant deviations from the SM predictions for gR and gi may be found. However, the 
SM values are only known for gR up to one loop in QCD while the prediction for gi is not 
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published. Moreover, most of the analysis frequently assume real values for and g^. The 
quantum corrections coming from the SM should be under known in order to discriminate 
the SM and new physics contributions from data. 

Let us briefly review the experimental status for the constraints on these tensor couplings. 
Indirect limits on gi and gR can be obtained from b ^ sj in the measured branching ratio 
B{B ^ Xsj). The results from a recent analysis [[T8l are given in the first line of Table 1. 

Table 1: Bounds on g^ and gi. The first line shows the indirect limits from b sy. The 
second and third lines are limits obtained from simulations for the LHC. The last two lines 
show 3 o discovery limits intervals: fourth line limits are from simulations for the LHC and 
the last one is from b sy. 



Reference 




gR bound 


gt bound 


in 


95%C.L. 


-0.15 <gi?< 0.57 


-0.0015 <gL< 0.0004 


in 


2o 


-0.026 <g/e< 0.031 


-0.058 <gL< 0.026 


m 


lo 


-0.012 <gR< 0.024 


-0.16 <gL< 0.16 




gR discovery limit 


gi discovery limit 


m 


3a 


\Re{gR)\> 0.056 


Re{gL) > 0.051 or 
Re{gL) < -0.083 


|/m(^^)| > 0.115 


> 0.065 


[fT3lfT8l 


3o 


Re{gR) > 0.76 
or 

Re{gR) < -0.33 


Reigt) > 0.0009 or 
Reigi) < -0.0019 


|/m(^L)| > 0.006 



The constraints on gt are much stronger than those on gR due to the chiral / nib enhance- 
ment factor which comes together with the gL coupling in the 5-meson decay amplitude. 
These bounds are obtained assuming that all anomalous couplings are real and that only one 
of them is non zero at a time. 

The top width Ft is an observable which is sensitive to the absolute strength of the tbW vertex 
but with no particular sensitivity to the anomalous couplings. Another test of the Lorentz 
structure of the tbW amplitude is the measurement of the polarized decay fractions Br{t 
bW.^) into bosons of helicity X = 0, =Fl-The SM values for this W-polarized widths 
are known up to one loop QCD and electroweak corrections [|8l and the contribution to the 
helicity fractions from the anomalous couplings defined in Eq.Hl were computed in ll9l[T0l. 
However, these fractions are sensitive only to ratios of the couplings. Other observables can 
be obtained from the single top production at the LHC [|6l. From a combined analysis of 
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the single top cross section and the three helicity fractions, the four anomalous couplings 
in Eq.[[Tl can be determined. Also, gR could be measured from the energy and angular 
distributions for polarized semi-leptonic and hadronic top-quark decays as was studied in 

m. 

The Tevatron found no deviation from the SM in top quark physics. It has also investigated 
the anomalous couplings based on the single top quark production cross section [|20l . These 
are the first direct experimental bounds but they are not competitive with the indirect ones 
already mentioned. In the near future the LHC high statistics data on top quark decays will 
allow the direct determination of the tensor couplings gi and g^ within a few percent level. 
In particular, in simulations for tt production and decay in dileptonic [[211 and lepton plus jets 
channels |[T2l[T9ll2TI . forward-backward asymmetries [|T2l and a double angular distribution 
in t-quark decay [|T9l were studied. In [[T2l . with only one non standard coupling different 
from zero at a time, intervals for detection or exclusion at two standard deviations (both 
statistical and systematic uncertainties included) for gi and gR were predicted for the future 
LHC data; they are shown in the second line of Table 1 . The LHC will possibly improve 
the sensitivity to by an order of magnitude when compared to the indirect bounds from 
sy. A combined fit, using the four couplings Vl, Vr, gR and gi as parameters, and taking 
into account the expected uncertainties at the LHC for top-quark decay was presented in [|6l 
and it is shown in the third line of Table 1. The sensitivity to gR shown in the second and 
third lines of Table 1 is similar to the results of [fT9ll2n where ti production and decay into 
leptonic and dileptonic decay channels were analyzed. 

Recently, new helicity fractions of the W where defined and investigated for polarized top 
decays. The spin matrix for polarized top decays was obtained in [[T3l : they also considered 
new observables derived from the normal and transverse polarization fractions. A similar 
approach in asymmetry observables was also studied in Tau physics in recent years [|22l . The 
three different sets of W helicity fractions defined for the polarized top quark were shown to 
open the possibility of new observables particularly sensitives to both the real and imaginary 
parts of the tensor couplings. They compute the 3a discovery limits for gR and gi assuming 
that they are either real or purely imaginary and allowing only one coupling to be different 
from zero at a time. The exclusion intervals are shown in Table 1 in the fourth line. As a 
reference for the comparison of the potential of the LHC they also derived the 3o discovery 
limits from b ^ syin lITSl ; this is shown in the last line of Table 1 . 

It is generally believed that beyond the SM theories will be probed in top quark physics. 
These theories, in general, not only will induce non zero values of gR and gi but also may 
be responsible for new exotic decay modes. The top dominant decay mode, t — )• bW^, was 
investigated in many extension of the SM. In particular, the decay rate and polarized decay 
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fractions were studied in many theories such as the Two Higgs Doublet Model (2HDM), 
the minimal supersymmetric SM (MSMM) and top-color assisted technicolor (TC2). These 
results were reviewed in [|71 where some of the anomalous couplings were computed for 
these models. For the first two of them they found the general feature that \gR\ ^ \gL\ and 
(sr) I ^ 1^"^ {sl)\- Besides, they found that values of gR up to 0.5 x 10^^ are possible for 
low tan P, while only 0.2 x 10"^ is expected for higher values of tan |3. Notice that these two 
last figures represent 8% and 3%, respectively, of the leading one loop gluon contribution. 
For TC2 models they showed that values for gR as big as 0.01 can be expected, and this 
represents 150% of the one loop gluon contribution. The general feature \gR\ ^ \gi\ and 
\R^{8r) I ^ {IinigL) I is also true for the SM, as it will be shown in the following sections. 



3 Electro weak corrections to the anomalous couplings gR 

and gL 

In the SM, at one loop, there is only one topology for the vertex correction diagrams that 
contribute to the anomalous gR and g^ couplings. This is shown in Figure 1(a) and we 
will denote this diagram as ABC using the name of the particles circulating in the loop. 
The QCD one loop gluon contribution (ABC = gtb) dominates the standard contribution in 
both the gR and gi dipole moments. All these diagrams, 19 as a whole, can be classified 
according to their dependence on the quark masses. As already mentioned, the tensorial 




(a) (b) 



Figure 1 : a. Topology of the one-loop SM Feynman diagrams for the quantum correction to 
the decay t — )■ bW^. b. Leading order diagrams for gR in the large limit. 

anomalous couplings we are interested in are chirality flipping magnitudes so, in general, a 
mass insertion is needed in order for the diagram to contribute. All contributions to gR need a 
nit mass insertion while the contributions to gi need a mj, mass insertion. The different mass 
insertions for each diagram is shown in Figure 2. Besides, some of the vertex have also a 
mass dependence. For all these diagrams there are three mass dependencies that are different 
for the case of gR and for gi. We use this fact in order to classify all the contributions 
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coming from all the diagrams. For there are two diagrams that have a leading mj-mass 

W W 




Figure 2: Mass insertions for the diagrams: in the left we can see the case for gj^ where a 
factor nit is present for each diagram while the ones at the right are for gi and proportional 
to m^. 

dependence, with non-decoupling in the ntj mass l|23l . They are the ones with thW and twoW 
circulating in the loop, where h is the Higgs boson and wq is the unphysical Z-boson; they 
are shown in Figure 1(b). These two diagrams have top mass-insertions that, together with 
the mass coming from the vertex, finally gives a mass dependence that is of the order 1 / r^, 
where = (mw/mt)^, with respect to the other diagrams. Next, there are 12 diagrams that 
also have a top mass insertion but do not have this 1 /r^ enhancement factor; all of them 
have a similar mass dependence. The QCD gluon exchange diagram also needs a top mass 
insertion for gj^. The remaining four diagrams do not have the 1 /r^^ enhancement factor 
but a suppression factor = {m},/m,Y coming from the mass dependence of the vertex 
and from a b-quark mass insertion. The particles circulating in the loop for these diagrams 
are bw^WQ, bw^h, htb and WQtb, where are the unphysical W-bosons. All these last 
diagrams turn out to be numerically insignificant as the small coefficient may suggest. 
Finally, the diagrams that dominate the final value for this electroweak correction are the 
ones of the first two classes we have already presented. In the computation we find that 
there is no numerical domination of the first ones with the mj non-decoupling effect over the 
ones without this effect. The gn coupling is finite and, up to one loop, the calculation needs 
no renormalization. However, not all the diagrams are infrared finite: some cancellation 
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occurs among some groups of them in such a way as to end up with a finite value. This 
fact was used as a check of our results. This is the case for the diagrams tyW and tWy and 
also for byW and bWy. There are no singularities when we sum each pair together for the 
computation. Some diagrams, like bWZ for example, contribute to the imaginary part of 
gR. In these cases we also used the Cutkosky rules and compared the result with the direct 
computation as a double check of our numerical calculation. Besides, some of the integrals 
can be done analytically up to the end and we verified the numerical evaluation of these 
results with the numerical evaluation of the Feynman integrals. This check was possible 
for the diagrams with a gluon or a photon circulating in the loop. All these facts allow for 
a multiple check of our calculations. As already anticipated, we can read in Table |2] that 
the value of the four diagrams with the factor are numerical negligible. Each diagram 



Table 2: Electroweak contributions to gR and gt. 



Diagram 


gR xl03 


gLXl03 


tZW 


-1.176 


-0.0141 


thW 


0.220 





tw w 


0.344 


0.0051 


thw^ 


0.462 


-0.0088 


tZw^ 


-0.050 


-0.0012 


tyw + tyw^ 


0.572 


-0.0094 


bWZ 


-0.623-0.664? 


-0.0201-0.0214? 


bWh 





0.0086-0.0120? 


bw^w^ 


(1.5 + 11.0?) X 10-4 


-0.0029-0.0167? 


bw^h 


(-4.3 + 8.6?) X 10-4 


-0.0019 + 0.0111? 


bw+Z 


-0.088-0.062? 


-0.00039-0.00028? 


bWy + bw+y 


0.114-0.509? 


-0.0270 + 0.0250? 


Ztb 


-0.397 


-0.0067 


ytb 


0.068 


0.0115 


Mptb 


-6.8 X 10-4 


-0.0109 


htb 


-6.2 X 10-4 


-0.0135 


T.{EW) 


-0.56- 1.23? 


-(0.092 + 0.014?) 


gtb 


-6.61 


-1.12 


Total 


-7.17-1.23? 


-1.212-0.014? 



contributes with a different sign to the final result so finally many of them are numerically 
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responsible for the final numerical value that can not be anticipating without an explicit and 
careful computation of all of them. The figures for each contribution of the diagrams to gR 
and gL is given in Tabled where we take the Higgs mass value m/, = 150 GeV. We can read 
in the Appendix A the expressions for each of the diagrams as well as the analytical results 
(Appendix B) for some of the diagrams for which it was feasible. We also show, in Appendix 
B, the limit m^, — )■ for these last formulas. This was also used as a check of our numerical 
evaluation of the integrals. The results are rather insensitive to the Higgs mass value in the 
experimental allowed Higgs mass interval [fT6l . The final result for the one loop electroweak 
correction for the magnetic right anomalous coupling is: 

gfM^ = -(0.56+ 1.23?) X 10^1 (2) 

Note that we have real and imaginary parts in this dipole moment and that the last is more 
than double of the first. These values are to be compared with the gluon contribution that is 
the dominant one: 

g| = -6.61 X 10-3. (3) 

This last result agrees with the one given in reference [ll51 if we put the numerical values for 
masses and couplings used at that time. The final result for the one loop computation in the 
SM is the sum of the last two values given in Eqs.jU and jS]]: 

= -(7.17 + 1.23/) X 10^1 (4) 

The real part for the one loop electroweak quantum correction for gR is 8% of the leading 
gluon-exchange contribution. The CP-even imaginary absorptive part, generated by elec- 
troweak corrections, may be measured with a similar set of observables as those considered 
in the literature to measure gR and, more specifically, Re{gR). Note that this imaginary part 
is 17% of the one loop Re{gf^) . 

The one loop electroweak correction for the gi coupling can be obtained in a similar way as 
in the previous calculation, however a b-quark mass insertion is present in all the diagrams 
because gt couples to a right b-quark. This factor dominates the numerical value of the final 
result for the gi electroweak contributions. As in the previous computation, there are also 
IR divergences in the same diagrams as in the preceding calculation and again they sum up 
to a finite result. The same checks we already explained before for gR have been used. We 
also find that an imaginary part for the electroweak contribution to gt shows up, so the final 
result is: 

gfw/ ^0.92 + 0.14/) X 10-4. (5) 
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This is to be compared and summed with the gluon contribution, that is real and the dominant 
one, in order to obtain the one loop result: 

^£ = -1.12x10-3. (6) 
The final result for the one loop computation in the SM is then: 

^f^ = -(1.21+0.01/) X 10^1 (7) 

We note that for gi the electroweak contribution is again 8% of the gluon contribution for 
gi, and the CP-even imaginary part has its origin in the electroweak diagrams. 



4 Conclusions 



We have computed the one loop electroweak values of the anomalous form factors gj^ and 
gi for the decay t — t- bW^. Both of them have real and imaginary parts. The imaginary 
parts come from the electroweak correction and, for gR, it is almost three times the real 
part, while for gi, it is 15% of the real one. Contrary to what happens in extended models, 
where the imaginary part are usually negligible if the new particles involved have higher 
mass scale than the top, we found that the absorptive parts of the dipole moments, which 
are induced in the SM by CP-invariant final-state re-scattering, has to be taken into account 
and may have physical effects that could be detected in the future through the observables 
proposed in the literature. Note that the gR coupling will be measured at the LHC and its 
absorptive contribution may be accessible in data and in the new observables defined in |fT3l , 
for example. For the SM one loop gR coupling, the imaginary part is about 17% of the real 
one while, for gi, it is only 1%. The value of the gi dipole moment, although proportional 
to m^, is only about one order of magnitude smaller than gR. The SM prediction for the 
real part of gL is Re^gi) — —0.0012 and is very close to the estimated 3o indirect discovery 
limit, Re{gL) < —0.0019, obtained from b sy in lfT3l based on the results of ifTSl . so 
any contribution coming from new physics that may show up may be in conflict with these 
bounds. Besides, the value of the electroweak corrections for these dipole moments, first 
published in this paper can, by themselves, explain deviations up to a few percent level in the 
observables, that are frequently discussed in the literature in connection to extended models. 
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A Expressions for the diagrams 



For the gR and gL couplings the contribution coming from each of the diagrams are written 
with the help of the following denominators: 

Az = xm{y-l)rl + l)y-rliy-\))-r^^ix-\) 

Bz=x{{ix{y-l) + l)rl+x-\)y-rl{y-l))-rlix-\)ixiy-l) + l) 
Cz = {x-l){xy~l)rl-rl{x-l)x{y-l)+r^^xy + xiy-l)ixy-l) 
{Ay,Bj,Cy} = {Az,Bz,Cz} (r, ^ 0) 
{Ah,Bh,Ch} = {Az,Bz,Cz} (r, ^ n,) 



The contribution of each diagram to is: 



tzw _^^Vbr-'fb f^^ ^ -2{at + vt)x{2{y-\)yx^-2yx^x+l) 



j dx j dy 
Jo Jo 



Az 



tjW _e^QtV;,r^,rb -2x {2{y - l)yx^ - 2yx + x + \) 



Sl 

jhW _ 

Sl - 

twow^ 

Sl 



sT' 



327^2 



X 



1287^2 



X 



dx \ dy 



/ dx \ dy 

Jo Jo 



-2x^y^ 
Az 



2 



128;r2„ e; 



128c2 7r2 



dx 



dx 







1 2x2 (43; -2) +2)3; 
^ Ah 
^ dy ^t + Vt){x- l)x 



Sl = X 



327^2 



dx 



Az 

'dy'-^ 



bWZ ^ ^^^tb fb 



dx I dy 



2x 



Qb {2{y-l)yx^ -5yx + x + A) + 



BzV 

Vb {2{y-l)yx^ + {y+l)x-2) 
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bWy_ e^QbVrlr..rh 2x {2{y - l)yx^ + {y + l)x - 2) 

~ 32712 Jo '^'^ Jo 

- 64nhl ^ Jo "^'Jo 

mn^r^sl Jo Jo Bz 
,t.n_ eXtrt 2x^{{x{y-l)-l)rl-x+l)y 

h^+z _ e^Vtlc^r^rb -2{ah -Vb)x^{y - 1) 

bw+y e^QbV;i,r,,rb 2x^{y-l) 
= 32^^/0 n'^^^ 

= i28ii;-4 '-^ ci 

while for we have: 

"128^^ 

/•I /•! 2x(vf (-23;;c2+;c+l) +af (-23;;c2 + 6y;c + ;c-5)) 

/ / J); — — 

h Jo Az 

tyw e^QtV;,r„ /•! /-i 2x{-2yx^+x+l) 
=^2^'' Jo''' Jo'' 

fdxtciy^-^^^ 
Jo Jo Ah 

jv Jo Ah 

('^KUw f\ . 2{at-vt){x-\)x 

^ X dx dy — -^^ — 

1284712 Jo Jo Az 

e^QtVrlr,. -2x{x-\) 
= -^2n^''Jo'''Jo''^^ ^'""^'^ 



thW _ ^^Kb 
tWQW _ __^__tb__ y 

1287r2r,,52 ^ 

thw- ^ ^ ^tb ^ 

1287r2r,,52 

tZw- 
Sr 



tyw 
Sr 
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bwz ^'^Vb^^ f\ /"^ 2{ab + Vb)x{2yx^ -{3y+\)x + l) 
=128^^/0 

Jo Jo 



- 32712 ^io^-^io"" By 

- 6471 V ""^-^ 



1287i2rH,4 Jq Jo ^ Bz 

- 12847^2 ^ io "^""yo 5^ 

fcw+Y_ e'6z.V;,rw /-l^ /-l^ -2x2(3;- 1) 
Sr — X 



/ dx j dy 
Jo Jo 



32712 -io"^7o"-' 5 



51247124 



y ^ y ^ -4(afc + Vfc)-^^ - 1 ) + 2) + v,x(3; - 1 ) )}; 

wofi ^^^^/fo'-^ /"^ , f^, -2rb{x-l)x{xy-l) 
12871-^^^, Jo Jo Cz 



Jo Jo 



1287r2rH,4 -/o Jo ""^ Ch 



B Some exact results 



The following integrals, corresponding to diagrams with a photon or gluon circulating in the 
loop, tyW, tyw^, bWy, bw^y, ytb and gtb, can be done analytically. Using the notation 



we have: 

jtyW+tyw' 



4 
A 

2 
rb 



8if 



'tb^w 



327^2 



Q X Itf 



1 + - 



Arl 



log 



i + ^f + A 



4rlA 



log 



A^ -A 



(l-'-^ + '-i + A)(l-r2+3r2 + A) , /I 



C + ^^ + A 

4 -4+ A 
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2/71 



[2 - 3rl + rl + rl + rl{\- IrD] - 2 + 2(2 - r^, + r^) log(2r2) + 



Arl 



jtWy+tw+y 



jltb 



^ 1^ 

with A = 
limit, as: 



jtjW+tyw 



yth 



jtWy+tw+j 
^R 



jtWy+tw^y 



jytb 



2rb 
A 



2 
A 

4rt 
A 



l-r^ + 3r^ + A 



(l-r2+r2 + A)2 
log(l 



A^ -A 



2 2 



+ A)] 



A^-A 



l-ri + ri + A 



^+rl-rl-A 
1 



log 



log 



A^ -A 
A^ -A 



1 — 2(r2 + r^) + (r^ — )2. These expressions can be written, in the — )■ 



l-r2 



+ 



3(1 



[3(r4 - 4r2 + 3) + 2(2 + 2r2 - r^) \og{rl)\ + 0(r4) 



-2/71 



2r 



[8r2-3r4-5 + 2(r2-2)log(r2)] + 
9r2 _ ,6 _ 8 ^ 3^^4 _ 2^2 _ 1) log(,2 )] + ^(,5^ 

.2 -| r 

l + (2-r2)log 



1 -r2 



0(r,^) 



[log(r2) - (3 - 2rl + r^) log(l - ) + (2 - 2r2 + /j log(r2 )] + 



4r? 



2rb 
1 -r2 



+ 



[log(r2) + (2 - rl) log(r2 ) - (3 - rl) log(l - ri 



2rl 



n [1 +21og(r2) +21og(r2) -41og(l - 4)] + 0(ri) 

-^log(l -4) + -T^^Tj [1 +log(r2) _21og(l -4)] + 0{rt) 



w 

1 -r2 
2r3 



log(r^ 



1+r: 



^log(l-ri 



(1-4)3 



[l+21og(r2)-41og(l-r2)] + 0(r5) 
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